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ABSTRACT 

We investigate the spatial extent and structure of the Pegasus dwarf irregular galaxy 
using deep, wide-field, multicolour CCD photometry from the Sloan Digital Sky Survey 
(SDSS) and new deep HI observations. We study an area of ^^0.6 square degrees centred 
on the Pegasus dwarf that was imaged by SDSS. Using effective filtering in colour- 
magnitude space we reduce the contamination by foreground Galactic field stars and 
increase significantly the contrast in the outer regions of the Pegasus dwarf. 

Our extended surface photometry, reaches down to a surface brightness magnitude 
' /Xr — 32 mag arcsec"^. It reveals a stellar body with a diameter of ~8 kpc that follows 

a Sersic surface brightness distribution law, which is composed of a significantly older 
stellar population than that observed in the ~2 kpc main body. The galaxy is at least 
five times more extended than listed in NED. The faint extensions of the galaxy are 
not equally distributed around its circumference; the north-west end is more jagged 
Kyi than the south-east end. We also identified a number of stellar concentrations, possibly 

, stellar associations, arranged in a ring around the main luminous body. 

■ New HI observations were collected at the Arecibo Observatory as part of the 

ALFALFA survey. They reveal an HI distribution somewhat elongated in RA and 
about 0?3 wide, with the region of highest column density coincident with the luminous 
galaxy. The HI rotation curve shows a solid-body rotation behaviour, with opposite 
ends differing by 15 km s^^. There is a stream to lower velocities about 5 arcmin from 
the centre of the galaxy. 

We were able to measure ugriz colours in a number of apertures using the SDSS 
data and compared these with predictions of evolutionary synthesis models. The results 
indicate that the outermost regions of PegDIG are 5-10 Gyr old, while the inner kpc 
contains stars ^1 Gyr old and younger. The colours correspond to K-stars; earlier 
subclasses are located in the innermost parts of the galaxy. PegDIG appears to be a 
relatively low-mass object, with a total dynamical mass of 3x10® Mq of which only 
30% in stars and 2% is in neutral gas. 

The extended stellar distribution, the appearance of faint light extensions, and the 
lack of low column density HI tails rule out a possible tidal origin or a ram pressure 
stripping scenario. We propose that PegDIG is a fairly recent acquisition by the Local 
Group, since it does not appear to be disturbed by interactions with other galaxies. 

Key words: galaxies: dwarf - galaxies: individual (Pegasus) ~ galaxies: structure - 
Local Group 
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1 INTRODUCTION 

The Local Group (LG), our immediate neighbourhood, is 
the group of galaxies in which our Milky Way galaxy formed 
and is evolving. The Milky Way and M31 are the two dom- 
inant spiral galaxies in the Local Group, and each is sur- 
rounded by an entourage of lower-mass companions. This 
kind of a bi nary structure is fou nd in nearby galaxy groups 
as well fe.g. iKarachentsev et al.i r2002a b. 2003). Apart from 
the late-type spiral M33, the more than fifty galaxies form- 
ing the LG are either dwarf elliptical (d E), dwarf sphe roidal 
(dSph) or dwarf irregular (dirr) galaxies. iGrebell (|200ll ) char- 
acterised galaxies as dwarf galaxies if their total B-band 
or V-band magnitude is fainter than —18 mag. While the 
primary distinguishing feature between dirrs and dSphs is 
whether they have or do not have HI, there are other distin- 
guishing features. dIrr galaxies typically have central surface 
brightnesses of ^v,o ^ 23 mag arcsec"^ and total Hi masses 
Mhi ^ 10^ M0. As their name suggests, their optical ap- 
pearance is irregular, which tends to be caused by scattered 
Hll regions for the more massive dIrrs. They are mainly 
found at larger distances from massive galaxies. dEs are 
spherical or elliptical in appearance, typically with Mv 
17 mag, and have ^v,o ^ 21 mag arcsec"^ and Mhi ^ 10** 
M0. Along with the dSphs, they are usually found in the 
vicinity of massive galaxies. dSphs tend to have Mv S5-14 
mag, ^ , v,o > 22 mag arcsec"^, and Mhi ^ 10^ Mq. See 
iGrelpell for further details. 

Dwarf galaxies are, in principle, simpler systems than 
large galaxies. However, those objects studied in detail 
showed that even dwarf systems are complex in terms of 
their exte nded star formation histor ies and chemical evolu- 
tion (e.g, lGrebel|[l997l : lMateolll99d ). Moreover, in the LG 
there is evidence for past interactions among some of its 
galaxies, including th e accretion of dwarf ga l axies by the 
two large spirals ( e.g, Ibata ot al. 1994, 2001; Yan nv et al 



20081 ) 



2005 ; IZucker et al.i,2004 : .Belokurov et al.. 2007a : BeU et al 



The LG member galaxies reside currently in a volume 
with a radius of ~ 1 - 1.2 Mpc and their integrated mass is 
esti mated to be between (1.3±0.3) x lO" to (2.3±0.6)xl0" 



M(T) ( Courteau fc van den Berri] Il999l : IKarachentsev et al.l 



l2002d ). The LG is the best-studied galaxy group thus 



far and includes t he fainte st dwarf g alaxies ever detected 
(e.g. lZucker et al.l i2006: B elokurov et aL, 2006. 2007E)- The 
dwarf, l ow-luminosity bodies in th e LG all contain old popu- 
lations l|Grebel fc Gallaghei|[200i) and may represent r elics 
from the re-ionisation epoch ( Gnedin fc Kravtsovll2006h . Fi- 
nally, there are persistent claims that dwarf galaxies are 
dark-matter (DM) dominated. Studying their properties in 
the LG, the place where the most detailed work can be 
done, is a step towar d understanding the nature of DM (e.g., 
iGilmore et al]|2007l ). 

It is important to determine how large a galaxy is in any 
observational band. The apparent size alters the baryonic 
mass estimate and shows how extended the stellar compo- 
nent is in comparison to the gas. This has also implications 
on the possible evolution of the object, because the larger 
a galaxy's cross-section, the higher are its chances of hav- 
ing interacted with other galaxies since its formation. The 
shape of the outer regions of a galaxy is a particularly sensi- 
tive indicator of past interactions; these tend to strip away 



stars and ISM creating tidal tails and spreading debris in 
the vicinity of a galaxy. However, it is very difflcult to de- 
termine the morphology of the outer extremely faint surface 
brightness regions for galaxies in the LG. Their closeness im- 
plies a large angular extent. If the coverage of faint surface 
brightness levels over large areas is required, the necessary 
observing time for direct, wide-field imaging becomes pro- 
hibitive. In addition, it is important to note that distorted 
galaxy shap es at low light levels can be an artefact of Pois- 
son scatter (|Martin. de Jong, fc Rixll2008l ). 

The availabilit y of imaging data from the Sloan Digita l 
Sky Survey (SDSS. lYork et al.Tl2000l : IStoughton et al.ll2002D 
offers a new and different way to approach this task. This 
data set is very uniform, covers a large fraction of the sky 
in five spectral bands, reaches objects as faint as 24 mag, 
and is publicly available. The SDSS imaging data have al- 
ready used to explore th e structural param e ters of other 
LG objects (e.g. , Drac o: lOdenkirchen et aL l_(|20C)lb'), And 
IX: IZucker et all (120041'); An d X: IZucker et al.1 (l2007i) ; Urs a 
Major: IWillman et al.l l |2005^ : Leo II: IColeman et al.l l|2007h . 
and many others). 

In this paper we present a new interpretation of observa- 
tional material of the Pegasus dwarf irregular galaxy (DDO 
216=UGC 12613, hereafter referred to as PegDIG) based on 
very deep and wide-field surface photometry using imaging 
data from the SDSS. PegDIG is one of the fainter LG mem- 
bers and contains very few stars younger than 100 Myr (see 
also iGaUagh er et al 1 19981 ). This object is not to be confused 
with the dwarf spheroidal galaxy Pegasus II = Peg dSph 
= Andromeda VI (Karachcntscv fc Kar achentseval Il999l : 
iGrebel fc Guhathaku rta 1999; lArmandroff et al.ll 19991 1 

We trace PegDIG across its entire angular extent in the 
optical to the lowest accessible surface brightness levels us- 
ing SDSS data, since all previous studies concentrate d either 
on st ars near the centre of the galaxy (e.g., Gallaghe r et al.l 
|l99^ or present ed surface bri ghtness profiles covering only 
small fields (e.g. . Ivan Zeel2000D . We combine these data with 
new H I maps derived from Arecibo Qbservator;J^ data c ol- 
lected for the ALFALFA survey (e.g.. iHavnes et al.|[200i ). 



THE PEGASUS DWARF IRREGULAR 
GALAXY 



As discussed bv iGaUagher et al.1 (|l998l) . PegDIG exhibits 
many of the properties typical of dIrr galaxies such as Hi, 
young stars, and Hll regions, and an irregular appearance 
due to recent star formation. On the other hand, it is at the 
faint end of the luminosity range of dIrrs, its outer isophotes 
are fairly smooth, its gas content is comparatively low, and 
its recent star format ion activity is low as well, prompting 
iGallagh er et al.l (|l998l ) to suggest that it may be a dIrr/dSph 
transition-type galaxy. 



^ The Arecibo Observatory is part of the National Astronomy 
and Ionosphere Center, which is operated by Cornell University 
under a cooperative agreement with the National Science Foun- 
dation. 
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Figure 1. The central part 13.5 X 13.5 arcmin) of combined 
g, r, and i image of Pegasus dwarf galaxy derived from a weigiited 
combination of individual single-colour images from tiie SDSS. 
North is up and East is to the left. The oblique faint line is a 
satellite track that was registered on the SDSS r image. At the 
adopted distance of 1.0 Mpc, l" = 4.85 pc. 



2.1 Optical data, morphology, star formation and 
distance estimates 

The basic parameters of PegDIG, compiled from the liter- 
ature, are listed in Table [T] There has been significant dis- 
agreement among the various distance estimators presented 
in the past, ranging from 1.75 Mpc derived from Cepheids 
to 760 kpc derived from the tip of the r ed giant branch 

i TRG B) in colour-magnitude diagrams. Ide Vaucouleurd 
197511 quoted a distance of only 170 kpc to PegDIG, con- 
sidering the galaxy one of our closest neighbours. The three 
most recent determinations all use ground-bas ed photomet- 
ric data and em ploy the TRG B (McC onnachi e et al.ll2005l : 
iTikhonov 2006a) or Cepheids (|Meschin et al.l 12009'') as dis- 
tance indicators. These studies converge to a value of about 
1 Mpc, s lightly larger than the distance of 760 ± 100 kpc de- 
rived bv lGallagher et al] (|l998l ) using data from the Hubble 
Space Telescope. These distances render PegDIG a periph- 
eral, potential member of the M31 subgroup within the LG. 
We adopt here a nominal d istance to PegDIG of 1 Mpc. 

iGallagh er et al.l l|l998l ) studied the central region of 
PegDIG using imaging in the B, V and I bands with the 
Wide Field and Planetary Camera (WFPC2) aboard the 
Hubble Space Telescope (HST). Their field encompasses re- 
gions of recent SF in PegDIG. These authors obtained the 
deepest photometry of PegDIG so far, reaching objects down 
to V~25. They identified a main sequence as well as blue 
loop stars younger than 0.5 Gyr. These young populations 
are clustered in two central clump s. Moreover, there are 
olde r, more widely distributed s tars. [Gallagher et al.l (|l998l ) 
and iDohm- Palmer et al.l (| 19981 ) concluded that the model 
best fitting their measurements has a relatively constant SF 



rate over several Gyrs, with a period of enhanced SF about 
2 Gyrs a go. The colours of the main-sequence stars mea- 
sured by [Gallagher et al.l (|l998l ) indicate a relatively high 
extinction of Ay ~0.47 or E{B — V)= 0.15, and the stellar 
models could match the observations only f or a distance to 
PegDIG of 760 kpc. [Gallagher et al.l (|l998l ) combined var- 
ious measurements including the Hi column density and 
IRAS dust maps in order to constrain the reddening and 
pointed out that these estimators do not yield consistent re- 
sults. The adopted extinction is one of the main reasons for 
thei r shorter distance. 

iKrienke fc Hodg3 (|200ll ) used the colours and red- 
shifts of backgrou nd gal a xies v isible in the HST/WFPC2 
data of iGallagher et al.l (|l998l ) to estimate the internal 
reddening in the central region of PegDIG. After correct- 
ing for the Galactic foregr ound reddening using the val- 
ues of ISchlegel et al.l ^9M), they found that PegDIG it- 
self contributes very little to the reddening and se e ms to 
have hardly any interstellar dust. iKrienke fc Hodg3 (|200ll ) 
derived an internal red dening of E(B — V) ^ 0.03 for 
PegDIG, con sistent with 'Gallagher ct alj (119981 ). The total 
reddening of Krienke fc Hodge (2001) is lower than that of 
IGallagher et al.1 (Il998 ') by ^ 0.05 in E(B - V). 

lYoung et al.l ( 20031 ) estimated a current star forma- 



tion rate of 3x10 ^ M© yr ^ from the Ha fiux , base d 
on the br oadband a nd H q imaging in Ivan "Ze3 (HqoO). 
iHunter fc~Elmegreenl (|2004) hsted PegDIG as having the 
lowest star formation rate among the 94 irregular galaxies 
of their sample (4.4 x 10~^ Mq yr~^) and the lowest star 
formation rate per surface area (8.3 x 10"'* Mq yr~^ kpc~^). 

NED lists a major axis of 5 arcmin (~1.5 kpc at 1 
Mpc) a nd a minor axis of 2.7 arcmin for PegDIG. Pfikhonovl 
(|2006al ) found a 5 kpc thick disk (diameter ~16 arcmin at 1 
Mpc) by tracing the red giant stars on images obtained with 
the 6-m BTA telescope and with HST. He wrote that the 
low-luminosity blue stars in PegDIG are scattered through- 
out the 5'x3'body of the galaxy. The HST image allowed 
the decomposition of the blue star distribution perpendicu- 
lar to the disk; these stars "virtually disappear at 200 pc" 
while the thic k disk can be traced vertic ally to ~1 kpc. 

Recently. iMcConnachie et al.l (|2007l ) compared the stel- 
lar structure of PegDIG derived from Johnson V and Gunn i 
images from the INT Wid e Field Camera with the H i maps 
produced by lYoung et al. I (|2003t ). They concluded that the 
Hi distribution is different from that of the "regular, el- 
liptical" stellar distribution, and interpreted that as strong 
evidence for ram-pressure stripping of the PegDIG by an 
intergalactic medium associated with the Local Group. 



2.2 The interstellar matter: HI distribution, mass, 
and dynamics 

iFisher fc Tuiwl (|l975l ) were the first to detect Hi emission 
from PegDIG and, based on its low radial velocity, suggested 
it is a Local Group member. The 21-cm line emission they 
measured had a width of 43±7 km sec~^. The neutral hy- 
drogen in PegDIG was mapped by Lo, Sargent, fc Young] 
(,.1993.) with the VLA, and by jlo fl^man ct al. (1996) with 
the Arecibo radio telescope. ILo et al.l (|l993h showed that 
the Hi is concentrate d in the main optical body and that 
it shows some clumps. iHoffman et al.l (|l996l ) found that the 
Hi distribution was asymmetric, with the Hi peak offset 
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Table 1. Basic Parameters of the Pegasus Dwarf Irregular Galaxy 



Parameter 




Value 






Reference 


(1) 




(2) 






(3) 


a (J2000.0) 




23''28™36.25= 








(5 (J2000.0) 




+14°44'34.5" 








Bt (mag) 




12.50±0.15 




van Zee 1 


'2000) 


(B - V) (mag) 




0.64±0.04 




van Zee 1 


'2000) 


(U - B) (mag) 




0.00±0.09 




van Zee 1 


'2000) 


as (arcsec) 




104.4 




van Zee 1 


'2000) 


Distance estimates 














Method 


Dista,iice M^odulus 


kpc 




Reference 




(1) 


(2) 


(3) 




(4) 




PN 




200-1600 


Jacobv & 


; Lesser (1981) 




TRGB 


Si26 


1700 


Hoessel ^ 


'z Mould (1982) 




TRGB 


Si25.5 


700-1200 


Christian & TuUy (1983) 




Cepheids 


26.1±0.3 


1700 


Hoessel ot al. (1990) 




TRGB 


24.9±0.1 


950±50 


Aoaricio 


(1994) 




TRGB 


25.13±0.11 


1060±50 


Lee (1995) 




TRGB, RC 


24.40 


760±100 


Gallagher et al. ( 1998) 




TRGB 


24.80±0.07 


919±30 


McConnachie et al. (2005) 




TRGB 


25.04 


1020±20 


Tikhonov (2006a) 




Cepheids 




1070±50 


Meschin et al. (2009) 



Note: The abbreviations used here are as follows: as is the exponential scale length of the disk fitted by Ivan Zed l|2000h 
to the B surface brightness profile; PN refers to a distance estimate using planetary nebulae; TRGB is a distance using 
the tip of the red giant branch in the colour-magnitude diagram; RC is the same, but using the "red clump" . 



by ~ 1' from the centre of the outer iso-intensity contour. 
The rotation curve derived from the Arecibo observations 
becomes flat at about ~4.5 arcmin SE of the kinematic cen- 
tre, which is itself offset by ^^2.6 arcmin from the optical 
centre. On the NW side the rotation curve falls off after 
peaking at ^^2 arcm in from the kinematic centre. Note that 
lYoung et al] (j2003l ) mention that the low velocity of PegDIG 
very close to zero - causes significant confusion with Galac- 
tic Hi, thus the interpretation of single-dish wide-beam Hi 
measurements of this galaxy may be problematic. 

iKarachentsev. Karachentseva. fc Huchtmeierl (|200ll ) 
observed PegDIG with the Effelsberg radio telescope as 
part of a program to map the local galaxy population. 
They give only an upper limit of 11 mjy from the Hi line 
flux, which they attribute to a low Hi content and not to 
confusion from local Hi. Note that the 21-cm half-power 
beam at Effelsberg is 9.3 arcmin. 

Ivan Zee! (|200Cll ) reported results of a study of isolated 
dwarf galaxies, where "isolated" implies a distance of 100- 
200 kpc from the nearest neighbour, using UBV and Hq 
imaging, combined with VLA H I mapping. PegDIG was in- 
cluded in her sample and yielded an integrated flux of 29.90 
Jy km s~^ and an Hi line width of 40 km s~^ for a recession 
velocity o f —183 km s~^, bas ed on the VLA observations re- 
ported in I Young ei~^ (|2003l ). van Zee's (2000) photometry 
showed that PegDIG has the reddest colours of all galaxies 
in her sample. Evolutionary synthesis models indicate that 
the stellar population could be the result of a single major 
star formation (SF) burst a few Gyr ago that ran out of 
material for further SF. 

Fur ther Hi mapp ing of PegDIG was done at Wester- 
bork bv lStil fc Israel (2002). as part of a project to study 
the hydrogen in a sample of dwarf galaxies. PegDIG was 



consider ed by these a uthor s as one of their most isolated 
objects. IStil fc Israeli (Hfloi') found a 9.0 mJy upper limit 
to the 1.4 GHz continuum flux from the galaxy in the 1.35 
MHz band, and a 21-cm line flux integral of 16.3±0.5 Jy 
km s~^ , implying a total H I content of 4x 10^ M0 at the as- 
sumed 1 Mpc distance. The map they presented in their Fig. 
4 shows the Hi concentrated at q(2000) = 23''26'"04' and 
(5(2000) = +14°27'"40", although the 13" synthesised beam 
was rather elongated in the north-south direction because of 
the declination of the object. The line profile shown in their 
Fig. 2 is single-peaked and very narrow, with a FWHM of 
less than 20 km s~^ and a weak ~ 80 km s~^ tail to higher 
recession velocities (the systemic Hi velocity of the galaxy 
being —189 km s"'^). 

The mo s t rece nt published H l map of PegDIG was by 
lYoung et al.l (120031 ') using the VLA and combining observa- 
tions performed in the C and D configurations. Their inte- 
grated 21-cm profile has a 50% full w idth of 24.6 km s~^ 
and a flux integral of 29.9 Jy km s~ M Young et all (|2003l ) 
display in their Fig. 5 a contour map of PegDIG with the 
lowest column density level at 10^^ atoms cm~^ and a ~10 
arcmin extent that shows the gas arranged in three main 
clumps. In addition, a region 1.5 arcmin northwest of the 
centre shows a double Hi proflle; this is interpreted as an 
expanding bubble with a radius of one arcmin (~200 pc at 
the 760 kpc assumed distance). 

'Jackson et al.l l|200(j ) presented spatially-resolved maps 
of PegDIG at 4.5 and 8 pm obtained with the Spitzer IRAC 
camera. These maps show only low surface brightness emis- 
sion from the galaxy, indicating low amounts of h o t dust 
grains and PAHs in the ISM. This, Ijackson et"al] (|2006l ) 
propose, is the result of the destruction of grains by super- 
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Table 2. Photometric parameters for the unresolved part of PegDIG in the SDSS 



Filter 




Circular Apertures 






2D GALFIT 






b/a 




a 


n 


b/a 




n 






[mag/sq.arcscc] 


[arcsec] 






[arcsec] 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


u 


0.40 


24.24±0.10 


89±9 


1.05±0.05 


0.40±0.01 


172: (98:) 


0.96±0.04 


a 


0.40 


23.10±0.05 


98±5 


1.03±0.02 


0.40±0.01 


172±1 (108) 


1.04±0.02 


r 


0.40 


22.73±0.04 


114±5 


1.10±0.02 


0.41±0.01 


172±1 (126) 


1.18±0.02 


i 


0.40 


22.46±0.05 


118±5 


l.lliO.03 


0.39±0.01 


177±1 (139) 


1.25±0.02 


z 


0.40 


22.40±0.07 


127±7 


1.33±0.06 


0.40±0.01 


164±2 (125) 


1.22±0.02 



Table 3. Total magnitudes of the Pegasus dwarf irregular galaxy 



Filter 


Circular Apertures 


2D GALFIT 






Apparent 


Total 


Total 






[mag] 


[mag] 


[mag] 


[mag] 


(1) 


(2) 


(3) 


(4) 


(5) 


u 


13.64±0.03 


13.60±0.46 


13.20±0.03 


0.34 


9 


12.26±0.02 


12.23±0.20 


12.01±0.02 


0.25 


r 


11.67±0.02 


11.63±0.17 


11.48±0.02 


0.18 


i 


11.35±0.02 


11.32±0.21 


11.16±0.02 


0.14 


2 


11.37±0.03 


11.32±0.29 


11.08±0.03 


0.10 


(« - 9) 




1.37±0.50 


1.19±0.04 








0.60±0.26 


0.53±0.03 




(r - i) 




0.31±0.27 


0.32±0.03 




(i^z) 




0.00±0.36 


0.08±0.04 




U 




12.88 


12.49±0.05 


0.36 


B 




12.68 


12.43±0.03 


0.28 


V 




11.87 


11.69±0.03 


0.22 


Rc 




11.41 


11.27±0.04 


0.18 


la 




10.89 


10.74±0.05 


0.13 



nova shocks and the inability of the ISM to regrow them in 
a regime of low or zero star formation rates. 



3 THE SDSS DATA FOR THE PEGASUS 
DWARF IRREGULAR GALAXY 

An image of PegDIG was created by extracting 
the relevant area from the SDSS data release DR5 
ijAdelman-McCarthv et al] |2007^ imaging data set. Owing 
to the distance of PegDIG, the high stellar densities in 
its central regions, and an average seeing of ~ 1.5" in 
the SDSS imaging data, individual stars can no longer 
be resolved as point sources in this area. We performed 
surface photometry separately for the inner galaxy parts, 
and for the outer parts, where stars still can be resolved. 
The two independently-derived surface photometry results 
were subsequently combined into a single, global surface 
photometry profile using areas of overlap. The results for 
the outer part were matched to those for the inner part, 
yielding a single surface brightness profile (SEP) per band 
for the entire galaxy. 



3.1 Integrated photometry and Surface 
Brightness Profiles 

PegDIG itself was not identified as a separate entity by the 
standard SDSS pipeline (see, for example, the SDSS DR6 
database). The integrated photometry, the creation of SBPs, 
and th eir analysis were done in the manner described in de- 
tail bv lKniazev et af] l|2004 ). First, g, r and i images for a 
wide region around the galaxy were extracted from the SDSS 
database. These were combined with weights to form an im- 
age of the object that is deeper than any of the single SDSS 
bands could offer and that has essentially the same stellar 
point spread function (FWHM~ 1") as the pipeline-reduced 
SDSS data. The central part of the combined gri image is 
shown in Figure [T] The galaxy location was defined above 
the 3cr noise level on the smoothed combined image, and all 
background sources were subtracted using SDSS database 
coordinates and using additional masking. The background 
was subtracted and SBPs for all ugriz filters were created 
in circular apertures with a uniform isophote step size of 2". 
The calc ulated profile s were fitted on a logarithmic scale, 
following [sirs3 l|l968l ): 



At(r) ^110 + 1.086 ■ (r/a)". 



(1) 



where R is the distance along the axis and n is an additional 
parameter. 

To check the stability of our results we performed ad- 
dition al two-dimension al modelling with the GALFIT pro- 
gram (|Peng et al.ll2002l ) using a one-component Sersic func- 
tion for all ugriz filters. To exclude all foreground stars from 
the fitting procedure we used the same mask image that was 
calculated on the previous step, where the method of circu- 
lar apertures was used. 

The results for both methods are summarised in Ta- 
bles [5] and [3] Table [2] presents all the model parameters 
for both methods, except for the total magnitudes. In case 
of circular apertures, the b/a axis ratio was calculated only 
once when building the mask for the galaxy location and 
was kept constant for the different filters. When GALFIT 
was used, the b/a axis ratio was allowed to vary as one of 
the parameters of the fitted model. GALFIT uses the fol- 
lowing form for the Sersic function: 



(2) 



where /Xe is the effective surface brightness and is the ef- 
fective (half-light) radius. Generally ct — Ve/bn, but b„ is dif- 
ferent for each n and can be only found numerically. A good 
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approximation that can be used i s b„ ~ 2/n — 0.324, but this 
is va lid only for 0.07 ^ n ^ 1 (|Truiillo. Graham, fc CaonI 
I2OOII ) and is not correct here, since n > 1. For this reason, 
our recalculations from to a are presented in brackets in 
column 7 (r^) of Table [1 The GALFIT solution for the u 
filter was unstable and for this reason we fixed Ve to the 
same value it was for g. This parameter is marked with a 
colon. 

Table |3] presents the total magnitudes. In column 2, the 
apparent magnitudes derived by integrating the luminosity 
within a circular aperture out to the limiting isophotal level 
of the SDSS data are listed. The total magnitudes presented 
in Table |3] were ca lculated extrapolatin g out to r = 00 using 
equation (9) from lKniazev et al.l (|2004l ) and the parameters 
of the Sersic function listed in Table [H The rather large er- 
rors of these magnitudes result from the big errors of the 
input parameters in Table [5] The total magnitudes from 
GALFIT are also a result of flux integration out to r = 00. 
In the case of GALFIT, the total magnitude is a free pa- 
rameter of the flt. We also presented UBVRcIc magnitudes 
that were recal culated from the uqriz m agnitudes using of 
equations from ljordi. Grebel. fc AmmonI (2006). The mag- 
nitudes in Table |3] are not corrected for foreground Milky 
Wa y extinction, b u t exti nction corrections calculated using 
the lSchlegel et al] (|l998l ) prescription are shown in column 
5 of this table. 

The photometric parameters calculated with either of 
the methods are very similar, but the total magnitudes re- 
sulting from GALFIT are systemat ically larger. C omparing 
with the integrated photometry of Ivan leiHooO) we note 
that our 2D GALFIT values for Bt, (B - V) and {U - B) 
are very similar to those of van Zee. 



3.2 The Resolved Outer Part 

3.2.1 Photometric selection of PegDIG stars 

Tracing out PegDIG to much fainter surface brightness lev- 
els than allowed by the direct surface photometry described 
in the previous section is based on selecting stellar objects 
that probably belong to the galaxy and on rejecting those 
that likely belong to the foreground. This method of em- 
pirical photome t ric fil tering was first described and used by 
iGrillmair et~al] lll995h. It was first implemented for SDSS 
data bv lOdenkirchen et al.l (l2001al) for the globular c lus- 
ter Pal 5 and was used by lOdenkirchen et al] (|2001lJ ) to 
study the Draco dwarf spheroidal galaxy. Variants of this 
photometric filtering of resolved stella r point sources in th e 
SDSS were also deve loped and used bvlRockosi et alj (|2002l ) , 
ISmolcic et al] l|2007l ). and lColeman et all l|2007fr 

First, we selected all point sources imaged by the survey 
in a ~0.6 square degree region centred on the Pegasus dwarf 
galaxy. Approximately 8000 point sources in this area were 
classified by the standard SDSS pipeline as stars. Since the 
SDSS detects only the most luminous stars in PegDIG, the 
number of identifiable stars belonging to the galaxy is quite 
limited. For this reason we did not use any additional selec- 
tion criteria, but rejected only the sources for which the pho- 
tometric error in any of the three most sensitive passbands 
(gri) was larger than 0.4 mag. The spatial distribution of all 
the selected sources has a mean stellar density of ~2.85 stars 
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Figure 2. i versus (g — r) colour-magnitude diagram for 
stars selected in an elliptical annulus with a semi-major axis 
of 0.05-0.15 deg, ellipticity 0.6 and PA=— 55°, centred on 
(a, (5)=(352. 15°, -1-14.743°). 



arcmin"^ for the foreground stars (all areas, excluding the 
central part of the field where PegDIG is obviously located) . 

The crucial point in constructing a high-contrast map 
for the outer part of the Pegasus system is to achieve an 
efficient discrimination between the tra cers of this system 
and the foreground field stars. Following lOdenkirchen et al.l 
l|2001bl ). we created a training set by selecting ~1500 stars 
from the central part of the PegDIG (defined as an el- 
lipse with a semi-major axis of 0°.15, ellipticity 0.6 and 
PA=-55°, centred on (a, (5) = (352.15° ,-^14.743°). 

In Fig. [2] we show a colour-magnitude diagram of the 
point sources that are candidate members of PegDIG as 
mentioned above. The diagram shows a fairly amorphous, 
extended plume of faint sources that resemble the distribu- 
tion of stars in o ther distant dirr galaxies from early CCD 
observations (e.g. lTosi et"al]|l99ll : lGreggio et al.lll993l ). The 
SDSS data share several characteristics with these early, 
ground-based CCD data: For a galaxy at the distance of 
PegDIG, they are close to the detection limit and they 
are strongly confusion and crowding limited. The extended 
plume of stars in the SDSS stellar point source catalog is 
composed of young main-sequence stars, blue loop stars, red 
supergiants, luminous asymptotic giant branch stars, and 
stars at the tip of the red giant branch. In other words, we 
are sampling the most luminous stars in PegDIG, which rep- 
resent a range of different stellar populations and ages cover- 
ing many billions of years. While these shallow data are not 
suited nor intended for detailed studies of the star formation 
history of PegDIG, they are very valuable for studies of its 
structure. The fact that stars on the upper red giant branch 
are included is particularly important, since this will be help- 
ful when trying to trace the faint outer limits of PegDIG, 
which are typically dominated by old(er) populations. 

Following t he procedure, described by 
lOdenkirchen et al.l l|2001al lbl). we calculated new colour 
indices ci and C2: 
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Figure 3. Separation of foreground Milky Way stars from the stars in PegDIG. Left: Density distribution fp of Pegasus stars in the (ci, i) 
colour-magnitude plane that is shown as contour plot. Stars were taken from an elliptical annulus as described for Fig. [2] Middle: Density 
distribution fp of field stars in the (cl,i) colour-magnitude plane. Stars were calculated outside of an ellipse with similar parameters 
and with a semi-major axis of 0.25 deg. Right: Lines of constant number ratio (the population contrast) s = fp/fp. 



Cl = 0.921 • 
C2 = -0.389 ■ 



'-r)+ 0.389 ■ (r - i) 
- r) + 0.921 ■ (r - i) 



(3) 



The re lations in the equat i ons [3] a re very close to those 
derived by lOdenkirchen et al] (I2OOIIJ) for the Dra c o dSph 
and similar to those derived by Odenkirchen et ahl (|2001al ) 
for Palomar 5. 

After constructing the colour index ci, which is a lin- 
ear combination of the {g — r) and (r — i) colours, the next 
step was to design an empirical filter mask in the (ci,i) 
colour-magnitude (CM) plane that will allow the optimal 
separation of the colour-magnitude distribution of sources 
belonging to the PegDIG population from that of the fore- 
ground field stars. The basic assumption is that the stars 
in PegDIG are of specific types, based on their colours and 
apparent magnitudes, and will occupy a specific locus in the 
CM diagram. This is shown in Figure |3] The region with 
maximal s values in the right panel of Figure |3]delineates the 
locus in the CM plane with the highest fraction of PegDIG 
stars relative to field stars. This last figure indicates that the 
stars selected from the SDSS data set as probable PegDIG 
members are mostly between 21sC i ^23.5 mag. 

The next step was to identify a level Sopt that yields the 
highest contrast for the selected area of PegDIG compared 
to the field stars. This implies that a filter mask for the 
selection of Pegasus stars should include as many points as 
possible with s ^ Sopt, where Sopt is a threshold value. To 
find the optimal number density threshold Sopt we computed 
the signal-to-noise ratio (S NR) for a range of s va lues, using 
the following equat i on from Grillmair et al.l (|l995h and from 
lOdenkirchen et al.l l|2001al lblh 



SNR(s) 



Np(s) - wNp{s) 
y/Np{s)+w^Np{sy 



(4) 



Here Np{s) is the total number of stars in the sample defined 
by s and for the region of PegDIG, and Np{s) is the number 
of stars in the same sample, but for the region where the 
foreground population is probed. The w parameter scales 
areas of these two regions. Since we are interested in the 
outer part of PegDIG, the SNR was optimised first for the 
elliptical annulus with a semi-major axis between 0.12 and 
0.20 deg from the centre of Pegasus. After that, other areas 



and different foreground regions were tested; we found that 
Sopt was practically identical in all those tests. 

The final selection filter constructed here removed 
~76% of the contaminating field stars and reduced the mean 
density of the contaminating foreground stars from 2.85 
stars arcmin"^ to 0.69±1.67 stars arcmin"^. The spatial dis- 
tribution of stars identified as foreground sources is shown in 
the top panel of Figure 3] It is clear that these stars are ran- 
domly distributed and do not show any strong concentration 
toward the bright core of PegDIG; a concentration could be 
expected if the filtering operation had been inefficient. The 
spatial distribution of stars with characteristics matching 
our final selection filter is shown in the bottom panel of Fig- 
ure |4l A visual comparison of the top and bottom panels 
of this figure shows that the objects selected by the filter 
as candidate PegDIG members do concentrate around the 
known galaxy, implying that the filtering operation indeed 
selected preferentially PegDIG stars very distant from the 
unresolved inner body. 

A stellar number density map of the central part of the 
studied field is shown in Figure [S] as isopleths. The surface 
density was derived through counts on a 30" by 30" grid and 
subsequent weighted averaging within a radius of one grid 
step. The thin lines show contours at the la level above the 
mean background density, where a is the rms of the back- 
ground stellar density fiuctuations, showing the lowest level 
at which the PegDIG stars start to be recognisable. Since 
any significant detection of the PegDIG population requires 
a surface density of at least 2a above the background, this 
level is marked by the thick contour in Figure [S] Levels of 
3(T, 5a and lOcr are plotted there with thick lines as well, 
to reveal the shape of the galaxy at different stellar num- 
ber densities; the plot shows that the overall distribution of 
PegDIG stars is approximately ellipsoidal, but that the outer 
contours appear deformed. It is also obvious that the NW 
end of the galaxy seems to be much more irregular than the 
SE one. These deviations, at the NW end, are very strong, 
showing up even at the 3cr contour. 



3.2.2 Structure of PegDIG 

To quantify the size, shape and orientation of the Pegasus 
dwarf irregular galaxy we fitted a two-dimensional model 
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Table 4. Parameter values for the best-fit model for the surface density distribution of PegDIG 
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Figure 4. The spatial distribution of SDSS selected point sources 
in the direction of the Pegasus dwarf galaxy. Top panel: Stars se- 
lected by the algorithm as belonging to the Galactic foreground 
appear randomly distributed in this image. Contours of equal stel- 
lar surface density at a level of Icr and 2a are shown with thin 
lines ((T is the rms background variation). There are some con- 
centrations at only a Icr level in the center. Bottom panel: Final 
spatial distribution of ~3400 photometrically selected stars in the 
PegDIG area. The stars selected as candidate PegDIG members 
are concentrated around the recognised location of the galaxy. 
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Figure 5. Top: Contour plots of the observed distribution of 
stars in PegDIG. Contours of equal stellar surface density are 
shown at 2a, 3a, 5a and lOcr with a thinner contour at la. The 
profiles in the centre of the galaxy show a "hole" produced by 
missing information in the high stellar density part of the galaxy. 
Dot-dashed lines show the contours of the best-fit exponential 
model (2D Sersic profile) at la, 3a, 5a and 10a levels. Note the 
"dendritic" extensions of the faint contours in the north-west part 
of the galaxy. Bottom: Residuals of the fit shown in the top panel, 
rescaled to the level of the mean background counts. As before, 
the contour with thinner lines corresponds to levels of ±la. Short- 
dashed (blue) lines show negative levels of the residuals. 
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Figure 6. Radial profile of tlie surface density of stars in PegDIG. 
Data points with error bars show the profile of the observed distri- 
bution (star counts in elliptical rings around the centre of PegDIG 
with parameters of the best-fit model shown in Table |4] and with 
the mean background 0.69 stars arcmin"^ subtracted). The red 
solid line shows the major axis model distribution for the GAL- 
FIT 2D Sersic model (see Table |4] for parameters). The blue long- 
dashed line shows the best-fit ID Sersic model, where the point at 
~12 arcmin was excluded since it represents the extensions visible 
in the bottom panel of Figure |5] When doing a fit with this point 
included, the result of the ID fitting is very close to the result 
of 2D with GALFIT. The green dotted line shows the major axis 
model distribution for the GALFIT 2D King model (see Table H] 
for parameters). Points located to the left of the vertical dotted 
line at R=4 arcmin were not included in the fit for either model. 
The vertical short-dashed line indicates a radial distance of ~450 
arcsec, where our SBP profiles in the ugriz filters terminate. 



to the observed surface density distribution. For that, the 
surface density was sampled on a 30" x 30" grid of non- 
overlapping cells. We performed two-dimensional modelling 
using GALFIT with a one-component Sersic function. The 
central part of the galaxy (which looks like a hole in our 
data) was excluded from the modelling with the special mask 
allowed by GALFIT. The best-fit parameters are given in 
Table U and the best fit is shown in Figures [S] and |B] The 
lower panel of Figure [5] shows the difference between the 
actual data and the fitted model. 

Although the statistical significance of any single fea- 
ture may not be extremely high, we note that there seems 
to be a ring-like distribution of peaks (clumps) around the 
unresolved part of the galaxy. The peaks in the ring configu- 
ration encircling the inner, unresolved part of the galaxy are 
highly significant, more than 5a above the background. The 
three peaks revealed to the north-west of the galaxy, between 
the la and 3a contours of the model, are of possibly even 
higher significance. These peaks co uld not be just artefacts 
of Po isson noise, as suggested by iMartin. de Jong, fc Rixl 
l|2008l) for other faint nearby galaxies, but here they are suf- 
ficiently significant, namely 3a or higher, as shown in the 
lower panel of Fig. [5] Such peaks are not seen at the oppo- 
site end of the galaxy and represent deviations of the faint 
outskirts of the galaxy from a pure Sersic profile, which are 
restricted to the NW side. 



As can be seen in Figure (5] the highest density of de- 
tected stars in the area of the peaks (clumps) is about 10-20 
a of the estimated background noise scatter. With the latter 
value 1.67 star arcmin"^ (see Section|3]2]), the 20 a level cor- 
responds to a star density of 33.5 star arcmin"^ or ~0.009 
star arcsec"^. In other words, the area with the highest den- 
sity of detected stars in the SDSS Pegasus data corresponds 
to about one star in a circle with a radius of ~6 arcsec. For 
an effective seeing of one arcsec and a pixel scale of 0.396 
arcsec there should not be any crowding problem in the re- 
gions with this and lower stellar density, either for detection, 
or for object classification as a star/nebula. 

To demonstrate the evidence for a much larger spatial 
extent of PegDIG indicated by our data, we show in Fig- 
ure|B]the observed radial stellar density profiles. These were 
derived from star counts in elliptical annuli, with parame- 
ters of the ellipses taken from Table [l] The logarithm of the 
mean density above background is plotted versus the radius 
of each annulus. The radius refers to the distance along the 
fitted major axis. The central 2.7' region was not included 
in the fit and is not shown. 

The best-fitting generalised Sersic model, described by 
the parameters listed in Table |4l is characterised by a radial 
scale length of a = 138 arcsec and an exponent n = 1.11. 
Comparing these numbers with those shown in Table [2] one 
can see that both the exponent n of the density distribu- 
tion and the scale parameter a are very close to the values 
derived from the gri images. As Figure [6] shows, the model 
profile fits our data very well over the entire region; there 
is no compelling evidence for a real cutoff in the radial pro- 
file of PegDIG within the current observational limits. Using 
SDSS data for the unresolved part of the galaxy we can trace 
the galaxy only up to a ~450" radius along the major axis. 
Using star counts, and improving the contrast of our density 
map, we can trace PegDIG up to ~800", but a distinct edge 
of the galaxy still remains undetected. 

The model fits provide a useful, smooth reference 
against which to identify peculiarities in the shape of the ob- 
served distribution of stars. Such peculiarities are the above- 
mentioned "ring" of stellar density peaks around the unre- 
solved part of the galaxy, and the three peaks at the NW 
end of the faint galaxy extension revealed by our filtering 
procedure and subtraction of the fitted model. We note that 
there are no counterparts to these peaks in the SE part of 
PegDIG. These features will be discussed below. 

3.2.3 Surface Brightness Profiles from Star Counts 

Our filtered star count data allow us to construct surface 
brightness profiles (SBPs). Such SBPs can be derived us- 
ing the filtered star count data, summing stellar magnitudes 
in elliptical annuli where the parameters of the ellipses are 
taken from Table |4j and normalising them to the area of the 
annuli. 

Assuming that the PegDIG stellar populations do not 
vary wildly in the outer parts of the galaxy (indeed they do 
not, as we show below), we can conclude that one sees the 
same profiles regardless of the method used. In the case of 
the unresolved part of PegDIG these were calibrated cor- 
rectly, but those constructed from the star counts were not, 
since stars below the detection limit of the SDSS would be 
excluded from the star counts. Using both data sets for the 
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Figure 7. Composite SDSS gr-band SBPs. The SBPs were calcu- 
lated in circular apertures for the unresolved part of the Pegasus 
dwarf irregular galaxy (blue points) and recalculated to the ma- 
jor axis using an cUipticity of e=0.61. The other data points for 
the SBP were calculated from star counts (black points) in el- 
lipses with parameters taken from the best-fit model for the stel- 
lar surface density map and corrected for the part of the colour- 
magnitude diagram that is not detected in the SDSS data. The 
red curves represent the best-fit Sersic model calculated from the 
SBPs for the unresolved inner part of Pegasus dwarf galaxy. The 
figures indicate that this model also fits the outer, fainter parts 
of this object reasonobly well. The SBP profile for the r filter was 
shifted by an additional —2 magnitudes for clarity. 
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Figure 8. Composite SDSS (g — r) and (r — i) colour distri- 
butions. One of the components of the colour distributions was 
calculated in circular apertures for the unresolved part of PegDIG 
(blue points) and recalculated to the major axis using an elliptic- 
ity e=0.61. The other component was calculated from star counts 
(black points) in ellipses with parameters taken from the best-fit 
to the stellar surface density map and corrected for the undetected 
part of the colour-magnitude diagram. The red curves represent 
the best-fit Sersic models that were calculated for the SBPs cre- 
ated for the unresolved part of Pegasus dwarf galaxy. The (r — i) 
colour distribution was shifted by -1-1.5 magnitudes for clarity. 



Dec, Index 59 




10 5 -5 

Hor, Po5, from Center (arcmin) 

Dec, Index 53 




10 5 -5 

Hor, Poo, from Center (aromin) 



Figure 9. Effective ALFALFA beams at two representative de- 
clinations. Top shows declination 14°47'"30'' and bottom shows 
14°41'"30''. The contours are spaced logarithmically, in steps of 
3.33 dB down from the beam maximum. 



part of the galaxy where the two methods overlap we can 
derive factors for the different filters and derive composite 
SBPs down to the faintest levels of PegDIG. These SBPs 
can be used to understand the evolutionary history of the 
galaxy, since the star formation in the outskirts of the galaxy 
presumably took place a long time ago. 

The interpretation of these correction factors is that 
they compensate the SBPs in the galaxy part created by 
using resolved PegDIG stars for the missing light produced 
by the fainter stars that are not detected on SDSS images. 
The underlying assumption is that the slope and cutoffs of 
the IMF are constant between the inner unresolved part of 
PegDIG and the outer part where the brighter stars are 
resolved and recorded. 

Composite SBPs in the gri bands were constructed from 
the unresolved part of the galaxy and from SBPs that were 
calculated from the star counts. Some of these composite 
SBPs are shown in Figure [7] Figure [8] shows the composite 
{g — r) and (r — i) colours. Both the [g — r) and (r — i) 
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Figure 10. Channel maps of the HI distribution and the integrated HI of PegDIG, derived from the ALFALFA observations. The 
central velocities of each channel map are (from left to right) —138, —153, —169 km s~^ in the top row, —184, —200, —215 km s~^ in 
the middle row, and —231 km s~^ in the bottom row. The total HI map is shown in the bottom-right panel, which also displays the 
beam used to produce all the maps. 



composite colours show very stable values for the outer parts 
of the galaxy, which justifies our previous statement. 

The correction factors for the gri filters were estimated 
visually, using both the SBPs and the colour diagrams. Their 
final accuracy is presumably about ±0.1 mag, since the 
colour diagrams are very sensitive to small photometric er- 
rors, and we estimate these to be ~0.05 mag for each band 
used here. The correction factors themselves are 2.25 mag 
for g, 1.85 mag for r and 1.65 mag for i filter. They were 
obtained by shifting vertically the values for the outer part 
of the galaxy to match those for the inner, unresolved part 
of PegDIG. 



range between -1600 km s"^ and 18,000 km The "fall" 
part of the survey maps the sky region from RA=22'' to 3*^ 
over most of the declination range accessible from Arecibo 
(5=0-36°). We collected the relevant HI observations cover- 
ing PegDIG from the ALFALFA archives as a grid contain- 
ing the galaxy. The ALFALFA observations are conducted 
at the Arecibo telescope with the 7-feed ALFA receiver. As 
described elsewhere (GiovaneUi et al. 2005a, 2005b, 2007), 
data acquisition for ALFALFA is done in a fixed azimuth, 
drift-scan mode. Each surveyed sky region is covered by 
two partly-overlapping drift scans collected at two differ- 
ent epochs in order to sample at better than Nyquist spatial 
frequency. 



4 HI OBSERVATIONS 

ALFAL FA, the Arecibo Legac y Fast ALFA extragalactic H I 
survey l|Giovanelli et al.ll2005^ is a blind neutral hydrogen 
survey that will ultimately cover 7074 square degrees of the 
high galactic latitude sky visible from Arecibo. It provides an 
extragalactic H I line spectral database covering the redshift 



The two scan passes are separated in time by several 
months, thus comparing the data from the two epochs helps 
in ruling out spurious detections. When data acquisition is 
completed over a given region of the sky, the individual drift 
scans are assembled and regridded to form three-dimensional 
data cubes or "grids". These grids are 2.4° x2.4° in size with 
1 arcmin sampling, and have 1024 channels along the spec- 
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Figure 11. Top panel: Position- velocity contour map along the 
major axis of PegDIG. Bottom panel: The global spectrum (to 
the 100 mjy km s~^ iso-intensity level) of PegDIG. The peaks 
between and -100 km sec"-*^ are Galactic hydrogen. 

tral axis. Since the velocity resolution is ~5.3 km s~^ (before 
Hanning smoothing), in order to cover the full ALFALFA 
velocity range which extends from -2000 to 18,000 km s^^ 
each grid is separated into four redshift subgrids, covering 
the velocity ranges: -2000 - 3300 km s~\ 2500 - 7900 km 
s-\ 7200 - 12,800 km s'^ , and 12,100 - 17,900 km 8"^ For 
mapping the HI in PegDIG we used only the first velocity 
grid of the data cube at the PegDIG coordinates. 

The central pixel of ALFA has low sidelobes, but the 
outer six pixels have significant coma lobes. Furthermore, 
the contributions of stray radiation and outer sidelobes to 
all the beams is not negligible. Data at each gridpoint in 
the ALFALFA grids are obtained by a weighted average of 
data from each beam that passes within a few arcmin from 
that point, thus the effective beam at any point has a quite 
complex shape. An IDL procedure for constructing effective 
beams for individual g ridpoints from the beammaps of the 
seven se parate beams Jlrwin et al. II2009I ) was written by J. 
Dowell (jGrossi et al. Ifiood ). Using this procedure, we have 
constructed effective beams at each declination pixel in the 
grid. Two of these are shown in Figure |9] 

The pronounced extensions shift in direction from one 
declination to another as little as 2' away as the dominant 
contribution to the received flux shifts from one ALFA pixel 
to the next. At a level of 10% of the main beam, these exten- 



sions to the effective beams cause strong, extended sources 
like PegDIG to have apparent emission reaching a few ar- 
cmin beyond the edge of the actual H I distribution at any 
declination where one of these extensions stretches toward 
the position of peak emission from the galaxy. 

Figure [TOl shows the individual HI channel maps (aver- 
ages of 3 ALFALFA velocity channels each) at the veloci- 
ties -137.9, -153.4, -168.8, -184.3, -199.7, -215.1, and 
-230.6 km s"\ The contours are drawn at 2.15, 4.64, 10.0, 
21.5, 46.4, 100, 215, 464, and 1000 mJy. The bottom-right 
panel in Figure [TO] presents the total HI map (integrated 
over velocity) with contour levels of 0.35, 0.74, 1.60, 3.5, 
7.4 and 16.0 xlO" atoms cm"^. The ALFA beam used to 
produce these maps is plotted at the lower-right corner of 
the total HI map; this is a fairly elliptical beam elongated 
north-south, with major axes of 3'.3x3'.9. The higher level 
contours agree well with the previous map of iHoffman et al] 
(| 19961 ) derived from Arecibo observations with the flat feed. 
The lower contours in our map show some differences, partly 
because the 1996 map was not a complete one with points 
spaced along major and minor axes as opposed to a full grid, 
and because the side-lobe structure is very different between 
the flat feed and the seven-feed ALFA. 

From an inspection of the gridview display of the indi- 
vidual channels in Figure [TO] we estimate that the structure 
in the lowest two contours to the north of the galaxy is 
possibly, as explained above, emission into the sidelobes of 
the beams that spanned that region and is probably spu- 
rious. The extensions to the east and southwest could also 
be sidelobe emission or could be real; this requires further 
observations. 

The global spectrum (to the 100 mJy km s^^ isophote 
level) for PegDIG, derived from the ALFALFA observations, 
is shown in the left panel of Figure 1111 The systemic HI 
velocity is —183.6 km s~^ and the width at 50% of the peak 
is 23.4±2.4 km s~^. The corresponding quantities at 20% of 
the peak are -185.2 and 38.6±2.4 km s"\ The HI profile 
is best-fitted by a Gaussian centred at —183.4 km s^^ with 
FWHM=23.6±0.3 km s'V The total HI flux is 28.1±0.1 
Jy km s~^, or 27.83±0.06 Jy km s~^ if the integral under 
the best-fit Gaussian is used. This corresponds to a total HI 
mass of 6.6T0^ Mq at the 1 Mpc nominal distance. There are 
no visible high- velocity HI wings. Note also that the signal 
from PegDIG is clearly distinct from the Galactic emission 
between —100 and km s~^ at this velocity resolution. 

Figure [TT] shows the major axis position- velocity (PV) 
contour map of PegDIG. To produce it we rotated the dec- 
lination (DEC) - right ascension (RA) maps by 25° (which 
we estimate to be the angle of the major axis of the HI dis- 
tribution away from the RA axis), then summed along the 
minor axis. The result is plotted as a contour map in the po- 
sition vs. velocity plane, with contour levels 21.544, 46.416, 
100.0, 215.44, 464.16, 1000.0, and 2154.4 mJy. The lowest 
two levels show the effects of the coma lobes, but the rest 
of the contours should be robust. The PV map indicates a 
solid-body rotation curve with opposite ends differing by 15 
km s~^, with a stream to more negative cz about 5 arcmin 
from the centre of the galaxy. The positive direction of the 
position axis is toward the SE. 

We measured a flux integral similar to that of Young et 
al. (2003) and almost twice that measured by Sill & Israel 
(2002). We conclude that, to a level of 2 mjy/beam, the 
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galaxy is 0.3 degrees wide in RA and 0.25 degrees in DEC. 
The peak of the HI is at RA = 23'' 28™ 55", DEC +14.75° in 
J2000 coordinates and at about -184 km s"\ The FWHM 
of the HI line is about 28 km s~^, about 10% wider that the 
width derived by Young et al. The HI distribution does not 
show an obvious compression ridge to the SE, or a tail to 
the NW smoothed out by our wide beam. The slight crowd- 
ing of the contours to the SE of the total HI distribution 
and slight elongation to the NW may be indicative of the 
compression ridge and tail inferred by Young et al. (2003) 
from their VLA map, but lacking a procedure to deconvolve 
the multiple asymmetric beams from the map we cannot be 
certain. However, we do not see any extension of the HI sig- 
nificantly beyond the outermost Young et al. contour even 
though ALFALFA is sensitive to column densities 4 or 5 
times lower than the sensitivity limit of the Young et al. 
map. 



5 THE STELLAR POPULATIONS 



Ivan Zed l|2000l ') presented B and {U - B), {B - V) profiles 
that extend only ~300 arcsec from the centre; these show 
that the stellar population changes from {B — V) ~0.6 one 
arcmin from the centre (/is —24) to ~0.8 some four arcmin 
from the centre. This is similar to the colour gradient found 
here for the unresolved part of the galaxy. 

Hierarchical galaxy formation scenarios predict that 
dwarf galaxies should show traces of an old stellar pop- 
ulation. And indeed, all nearby dwarf galaxies show 
clear evidence for the pre sence of old stellar populations 
l|Grebel fc Gallagher! 120041 '). i.e., populations older than 10 
Gyr, as traced by, e.g., the presence of a horizontal branch. 
This may also be the case for PegDIG, where traces of 
present-day very low le vel star formation (4 . 4xl0 ~^ Mq 
yr"""^) were found (e.g., iHunter fc ElmegreenI |2004| ) while 
most of the stars belong to earlier generations. Indeed, the 
deepest existing HST data show clear evidence for substan- 
tial populations with ages of 9 Gyr and younger as trace d 
by a very prominent red clump (Gallagher et al.| Il998h . 
Those data approach the detection limit just below the 
intermediate-age red clump, making it very difflcult to es- 
tablish the presence of a horizontal branch (although there 
are indications of a density enhancement in the correspond- 
ing region in the col our-magnitude diagram in Fig. 15 of 
ICallagher efaLlliggsl ). 

van Zed (|2000h found that PegDIG had the reddest 



{U — B) and (B — V) colours of all the galaxies in her sam- 
ple. She could not fit these with a constant star formation 
model, but could fit them with a continuous but declining 
star formation rate, or with an aging, a few Gyr old star- 
burst. She proposed that the very low star formation rate, 
now relegated to the innermost parts of the galaxy, might re- 
sult fr om the depletion of the galaxy's HI reservoir. Ivan 

mentioned also that the HII regions detected in her 
Ha images are so faint that they would not have been de- 
tected if PegDIG were more distant. 

To estimate the distribution of stellar population ages in 
PegDIG, we compared the derived integrated colours in its 
differ ent parts with model track s from the PEGASE2 pack- 
age (|Fioc fc Rocca-Volmeran"g3 1 19991 ) for metallicity val- 
ues of 2=0.02 (solar), z=0.008 (1/3 of solar) and 2=0.004 



(1/5 of solar). Since the photometric systems {u ,g' ,r' ,i' ,z') 
used for the calculations of PEGASE2 evolutionary tracks 
and {u, g,r,i, z) used in the real SDSS observations are 
sligh tly different, we ap plied the transformation formulae 
from iTucker et al.l (I2OO6I ) in order to correct theoretical val- 
ues to the {u, g,r,i, z) system. PegDIG has been found to 
have a low present-day ISM metallicit y of ~l/5 of solar 
(Skil lman. Bomans fc Kobuln ickv 1997 ) and poss ibly has 
more metal-poor older stars (Gallagher et al.lll998l l. Based 
on the models we used, tracks of even lower metallicities 
(<l/5 of solar) would not match the red colours observed 
in the outer parts of PegDIG and thus these tracks are not 
plotted. 

In Fig.[T5]we plot model tracks of colour evolution in the 
(g^r)o versus {u—g)o and (r — i)o versus {g — r)o diagram for 
a standard Salpeter IMF, with lower and upper mass cutoff 
limits of 0.1 and 120 Mq. All observed colours are corrected 
for Milky Way foreground extinction. Different tracks in Fig. 
1121 represent the colour evolution for continuous SF with 
constant SFR (dashed lines) and for an aging instantaneous 
SF episode (solid lines) as two extremes of all possible SF 
histories. The hexagons on the evolutionary tracks with the 
respective numbers mark ages in Gyr since the beginning of 
SF. 

We seek to constrain the mean age of the dominant 
population in different annuli by trying to fit simultaneously 
both colour-colour plots, the (u — g)o vs. {g — r)o and the 
(g — r)o vs. (r — i)o colours. Clearly, younger populations 
are more prominent in the central regions, where the mean 
colours agree relatively well with the 1 Gyr tracks. In the 
outermost regions the averaged age is 5-10 Gyr or less. Our 
data do not allow us to prove or to disprove the presence of 
an old population (older than 10 Gyr), but the comparison 
of the integrated colours with population synthesis models 
suggests that the contribution of such an old population to 
the integrated light is fairly minor. Hence the colour gradient 
appears to be consistent with an age gradient in the sense 
that younger populations are more centrally concentrated, 
while older populations show a more extended distribution. 
Of course, independent of the integrated colours, we know 
already from published earlier studies that an age gradient 
is present, as indicated by, e.g., the centralised distribution 
of the HII regions. 

As is well-known, colour gradients may b e caused by ei- 
ther a ge, metallicity, or reddening gradients ([Harbeck et al.l 
I2OOII : iHunter fc ElmegreenI l200q ) . We cannot unambigu- 
ously disentangle these three effects in our data, but we 
can at least attempt to qualitatively comment on their im- 
portanc e. Since the intrinsic red dening of PegDIG seems to 
be low (|Krienke fc Hodgdl200lh , we discard extinction as 
a significant contributor to the observed colour gradients. 
We have insufficient data to constrain a potential metallic- 
ity gradient, but we note that the slope of the metallicity 
gradient has to have an opposite sign to the observable one: 
the more metal-poor track has bluer colours compared to 
the more metal-rich track. In other dwarf galaxies for which 
more detailed data are available, metallicity gradients tend 
to be small despite considerable scatter at a given age (e.g, 
iKniazev et alll2005l : [ciatt et al.ll2008l : lKoch et al.ll2006l ). Al- 
together, it seems likely that our large-scale colour gradients 
are indeed primarily driven by age, although an age gradient 
may certainly also be linked with a metallicity gradient. 
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Figure 12. Two-colour diagrams with theoretical tracks from PEGASE2 for evolving stellar populations. The observed colours are 
corrected for Milky Way foreground extinction and are shown as red filled circles with error bars for the following major axis distances: 
0-50, 50-100, 100-150, 150-250, 250-400, 400-600 and 600-800 arcsec. {u - g)o could be only calculated up to a major axis distance 
of 400 arcsec. The observed colours were averaged inside each major axis distance bin and the ±lcr errors also reflect the scatter of 
colour in each distance range. Redder colours correspond to an increased major axis distance: the lowest-leftmost point represents the 
centre of PegDIG and highest-rightmost one corresponds to the outermost part. Top panels: Tracks with a standard Salpeter IMF, for 
instantaneous (solid) and continuous (dashed line) SF laws using different metallicities: blue - 2=0.004 (1/5 Zq), green - 2=0.008 (1/3 
Zq) and black - Zq. Filled black hexagons along the tracks, with respective numbers, mark the time since the beginning of SF: 1, 5 and 
10 Gyr. Bottom panels: Tracks for instantaneous (solid) and continuous (dashed line) SF laws with a standard Salpeter IMF (blue) for 
2=0.004 and non-standard Salpeter IMFs: slope is —1.6 for green and —2.85 for black. 



Population gradients have been id entified in early-type 
and late-type dwa r f galaxies alike (e.g. [ Harbeck et al.ll200ll: 
Parodi et al]l2 002': 'Makarova et al.l l2002l . l2005l : iLisker et al.1 
20061 : iHunter fc Elmcgrccn 200^ although not all dwarf 
galaxies show such gradients. Generally seen as colour gra- 
dients, such variations are usually interpreted as age gradi- 
ents. They appear to indicate longer-lasting star formation 
in the deeper, inner parts of the potential well where the 
star- forming material accumulates more easily and/or can 
be retained more easily. PegDIG seems to fit this trend as 
well. 



6 DISCUSSION 

One of the results of the present work is the deriva- 
tion of surface brightness profiles for PegDIG that ex- 
tend about 1.6 times further out than what was foun d 
in previous surface photometry studies (jTikhonovll2006al ). 
Other deep optical surveys of dirrs in the Local Group 



have found additional evidence for stellar distributions much 
more extended than pr eviously thought: e.g., NGC682 2 
(|de Blok fc Wa hei* '200d ), Leo A ('Vansevicius et al.' '200 ' 
and IC 1613 (Hattinchi e taLll2007l ) . Tikhonov (2005, 2006 
also found such an effect for many nearby face-on and edge- 
on dIrrs. 

Our deep surface photometry, using stars identified as 
belonging to PegDIG and eliminating Milky Way stars in the 
foreground, showed that it is possible to fit a Sersic surface 
brightness distribution to the galaxy at least to a radius 
of 14 arcmin~ 4 kpc. The stellar distribution shows some 
concentrations around the unresolved 2 kpc part, and some 
irregular extensions at the NW end of the major axis. Note 
that this region shows less HI than the SE part, when our 
maps are compare d with those for neutral hydrogen from 
lYoung et al.l ()2003l ). The HI column density shown in their 
Fig. 5 has a peak of ~ 10^^ atoms cm'^ at RA = 23''28™35", 

DEC= -f 14°44'"15", identical to the peak of the ALFALFA 
maps. 
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Our results allow the derivation of some general prop- 
erties of PegDIG. The total HI flux and the total B mag- 
nitude yield M(HI)/Lb ~0.4. A rough estimate of the to- 
tal dynamical mass can be derived from the PV plot in 
the left panel of Figure 1111 Assuming that the outermost 
gaseous regions of PegDIG rotate with ~60 km s~^ at a 
galactocentric distance of 4 kpc, the total dynamical mass 
is yidy-n — 3.3 X 10* M0. The neutral hydrogen contributes 
only 2% of this mass. The mass in stars can be roughly es- 
timated from the SDSS colours plotted in Figure [T2l these 
correspond to main-sequence K stars. If the light were pro- 
duced solely by K5 stars, the total mass in stars would be 
M, = 1 X 10* M0, about 30% of the total dynamical mass. 
This indicates that while PegDIG has a signiflcant amount 
of dark matter, it is not a dark-matter-dominat ed galaxy, 
con trary to the asser t ion of lAparicio et al] |l993). 

I Wilkinson et al.| (|2004l ) interpreted the sharp drop in 
the surface brightness proflle of the Draco dSph as the sig- 
nature of a kinematically cold stellar population in the outer 
parts of the galaxy. They proposed two possible explana- 
tions for this phenomenon, which they presented along with 
arguments against accepting them. One is a two-population 
model with a hot bulge and a cold disk or halo that does not 
fit either the observed kinematics nor the light distribution. 
The other is tidal reshaping by the Milky Way galaxy, which 
requires very tight constraints on the orbit of the dwarf 
galaxy. We mention this here since a similar phenomenon 
may be present in PegDIG, even though PegDIG is more 
than five times a s distant from the n earest massive spiral as 
the Draco dSph jGrebel et alll2003l ). 

iKarachentsev et al. I Hooi) mentioned that the major 
disturber galaxy for PegDIG, in terms of tidal interactions, 
is M31, with a 6 parameter of 1.2. This parameter is the 
tidal index, which is defined as: 



+ C 



(5) 



Here i is the index of the specific galaxy 6 is calculated for, 
k is the index of any other galaxy, is the mass of the 
fc-th galaxy and Dife is th e 3D space distance between th e 
z-th and the fc-th galaxy (jKarachentsev fc Makarov|[l999l l. 
Galaxies with 6 ^0 can be considered to be undisturbed, 
while objects with 6 >5 are considered to be highly dis- 
turbed. Based on the calculated value, PegDIG could be 
somewhat affected by M31, but not b y an extreme tidal 
interaction. iHunter fc ElmegreenI (|2004l ') listed M31 as the 
nearest neighbour to PegDIG, at a projected distance of 
450 kpc and a relative recession velocity of -1-117 km sec~^. 
Note though that since the actual orbit of PegDIG is not 
known, and it could even be on a fairly radial orbit around 
M31, a relatively strong tidal interaction in the past cannot 

be e xcluded. 

iMcConnachie et all (|2007l ') plotted some of the PegDIG 
neighbours in their Fig. 1 (lower right panel). Their argu- 
ment is that a comparison of the HI contours with the stellar 
distribution fits a morphology of gas being stripped away by 
ram pressure. They identify a "compression front" on the 
SE end of the galaxy. In order for ram pressure stripping 
to take place, they require the Local Group to be filled by 
a tenuous (~ 10~^ cm~^), rather warm (~ lO^K) gaseous 
medium. Ram pressure stripping would also modify the gas 
distribution in and around the galaxy. In fact, this is the 



main argument used by McConnachie et al. in bringing up 
the stripping possibility. The morphological modifications 
that should appear if the ram pressure stripping argument 
is valid are a leading edge compression region and the cre- 
ation of a tail of stripped material following the galaxy. 

The deep HI maps shown in Figure [10] lack the signa- 
tures of either a strong tidal interaction or of ram pressure 
stripping of the gas from the galaxy. On the contrary, the 
low column density gas observed at the outskirts of PegDIG 
and the lack of disturbances in the distribution of this tenu- 
ous gas imply that any presumed interaction with an inter- 
galactic medium, as propose d by McConnachie e t al (2007), 
can probably be discounted. iGrebel et al.l ()2003l ) presented 
simple estimates of the efficiency of ram pressure stripping 
by a homogeneous Local Group intergalactic medium and 
concluded that the densities are too low to have a notice- 
able effect. However. iGrebel et al.l (|2003l ) also argued that a 
putative clumpy medium could be rather effective. But the 
absence of disturbances in the PegDIG HI contours does not 
support stripping by a clumpy intergalactic medium. 

With the pre sent star formation rate derived by 
lYoung et al.l (Hooi) and the total HI mass measured in 
this work, the star formation could last for another 2x10^^ 
years. Such slow, continuous star formation extending over 
a Hubble time or more is typical for dirr galaxies (e.g, 
iHunter fc GaUagheJ [l98^ : [raxi .Zoc 2001,). The central HI 
column density is a bit short of the canonical threshold for 
star formation, but this ca nonical limit is not always valid 
for low-mass galaxies fe.g.. [Hunter et al]|l998l ). Is it possi- 
ble that we are seeing PegDIG during an extended lull in 
star formation, while the gas is settling back after having 
been distended by, e.g., supernovae in the last star forma- 
tion episode, akin to scenarios described by, e.g.. lDong et al.l 
(2003)? Amplitude variations in the intensity of star forma- 
tion are commonly observed in dIrr galaxies, so here also 
PegDIG wou ld conform to the typical pr operties of these 
objects fe.g.. iTosi et"allll99ll : iGrebelll 19971 ). 



iBegum et all l|2006l ) discussed the SF threshold in very 
faint low-mass galaxies based on synthesis observations with 
the Giant Metrewave Radio Telescope (GMRT). They found 
that while current star formation (as traced by Ha emission) 
is confined to regions with relatively large HI column den- 
sity [Nh/ >(0.4-1.7)xl0^^ cm"^], the morphology of the Ha 
emission is generally not correlated with that of the high HI 
column density gas. A high gas column density may be a nec- 
essary condition for star formation, but it is not sufficient, 
for their sample at least, to ensure that star formation does 
in fact occur. 

We can also rule out tidal deformation, since such dis- 
tortions are expected to be symmetric with respect to the 
centre, while for PegDIG any distorsions are relegated to 
the NW side. Our findings imply that PegDIG might be a 
recent acquisition of the Local Group, now in the outskirts 
of the LG and far away from any nearby massive galaxies, 
that was formed in a comparatively empty region without 
major disturbers. However, without knowledge of its orbit 
we cannot rule out past interactions with M31. 
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7 CONCLUSIONS 

We analysed images of the Pegasus dwarf irregular galaxy 
collected by the SDSS survey and found that the unresolved 
part can be fitted by a Sersic intensity profile down to ~ 30 
mag arcsec"^. Using very effective filtering in the colour- 
magnitude space of SDSS data, we reduced the contamina- 
tion by foreground Galactic field stars and significantly in- 
creased the contrast for the outer part of the Pegasus dwarf 
where we identified resolved stars that belong to PegDlG. 
This allowed the extension of the surface photometry to 
much fainter levels. 

Our extended surface photometry, reaching down to a 
surface brightness of /ig ~33 mag arcsec"^, revealed a ~8 
kpc wide stellar distribution following the same Sersic pro- 
file as found for the unresolved part, composed of a stellar 
population similar to that in the ~2 kpc main body, but sig- 
nificantly older. The distribution of colours across the galaxy 
body shows that the innermost parts of the galaxy contain 
the youngest population. A comparison to population syn- 
thesis models suggests a mean age of ~ 1 Gyr for this part, 
and we know from earlier work at HST resolution that the 
youngest stars have ages of only a few 10 Myr. The outer- 
most parts of PegDIG are much older with a mean age from 
integrated colours of at least 5 Gyr, indicating that there the 
contribution of younger populations is comparatively small. 
The total dynamical mass of the galaxy is ~ 3 x 10* Mq , of 
which about 30% is in stars and only 2% is in HI. 

We found that the stellar distribution of PegDIG is con- 
siderably more extended than previously thought. Mapping 
the HI distribution to very low column density levels at 
Arecibo revealed that the hydrogen distribution is slightly 
smaller than that of the stars revealed by the present study. 
PegDIG is therefore yet another dirr where the HI is not 
much more extended than the stellar distribution, as it used 
to be in the classical picture of dirrs. 

Our deep HI map shows an extended and fairly regular 
gas distribution with solid-body-like rotation. We do not 
observe low column density tails extending beyond the edges 
of the neutral gas distribution shown in earlier synthesis 
array images (Young et al. 2003). Tidal stripping therefore 
seems unlikely. On the basis of the HI observations alone, we 
cannot rule out ram pressure interaction with extragalactic 
gas, and the relatively small extent of the HI vis-a-vis the 
distribution of stars may even support the hypothesis that 
the outermost gas has been stripped. 

We identified faint extensions of the optical light distri- 
bution of PegDIG at the north-west end that do not follow 
the expected distortions caused by a tidal interaction. We 
also showed that a number of stellar concentrations - pos- 
sibly extended stellar associations - are located around the 
unresolved central galaxy body. Rings of stellar associations 
have been found in a number of dIrrs and could be a possi- 
ble sign of outward-propagating star formation but, on the 
other hand, PegDIG also has very young s tars in its inner- 
most regions as shown by the HST data ([Gallagher et al.l 
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